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We study a spin-1/2 model with triangular A'XZ-clusters on the orthogonal-dimer chain in the 
presence of an external magnetic field. First, we discuss the case where the triangular clusters are 
coupled via intermediate "classical" Ising spins. Diagonalization of the triangular XXZ-clusters 
yields the exact ground states; finite-temperature properties are computed exactly by an additional 
transfer-matrix step. A detailed analysis reveals a large variety of ground states at magnetization 
M equal to fractions 0, 1/4, and 1/2 of the saturation magnetization M = 1. Some of these ground 
states break translational symmetry spontaneously and give rise to doubling of the unit cell. In a 
second part we present complementary numerical data for the spin-1/2 Heisenberg model on the 
orthogonal-dimer chain. We analyze several examples of T = magnetization curves, entropy as a 
function of temperature T and magnetic field, and the associated magnetic cooling rate. Compari- 
son of the two models shows that in certain situations the simplified exactly solvable model yields 
a qualitatively or sometimes even quantitatively accurate description of the more challenging quan- 
tum model, including a case which may be relevant to experimental observations of an enhanced 
magnetocaloric effect in the two-dimensional compound SrCu2 (603)2. 

PACS numbers: 75.10.Pq; 75.30.Kz; 75.40.Cx; 75.40.Mg 



I. INTRODUCTION 



There is a two-fold interest in the statistical mechanics 
of many-body lattice spin systems. On the one hand, one- 
dimensional integrable many-body models"'^'^ have been 
intensively investigated, if not since the seminal 1931 pa- 
per of Hans Bethe,'^ then at least since the 1960-s. A 
linear chain of localized spins vifhich interact with each 
other via exchange interactions (Heisenberg chain) is the 
key component of almost all such models. The research 
on integrable spin chains is often considered to be sepa- 
rated from the underlying physical context, as their inte- 
grability is deeply related to modern algebraic structures 
(see for instance Refs. 4 and 5) and attracts mainly math- 
ematical physicists. On the other hand, various lattice 
spin systems serve as the basis for physical models for de- 
scribing magnetism of materials and related phenomena. 
As a rule, quantum lattice spin models corresponding 
to real materials rarely admit exact solutions. There- 
fore, one often has to resort to numerical calculations 
in order to obtain reliable results for the thermodynamic 
properties of such strongly correlated magnetic materials. 
However, for one-dimensional integrable models, there 
are powerful mathematical tools allowing one, at least 
in principle, to compute the thermodynamic functions 
exactly.^'^ Despite their impressive efficiency, these meth- 
ods are applicable only if the underlying lattice model is 
integrable, which is not always the case for real magnetic 
materials. However, there is another class of applica- 
tion for so-called classical one-dimensional lattice mod- 
els, namely as effective descriptions of the conformations 
of macromolecules.^ In that case, one can obtain an an- 
alytical expression for the partition function and, thus, 



for all thermodynamic functions. 

Lattice spin models with quantum and classical 
spins, ^ respectively Ising and Heisenberg^^^^ or Hubbard 
interactions^°'^^ combine properties of both quantum and 
classical spin systems. They can be solved exactly by 
the classical transfer-matrix method, while the entries 
of the transfer matrix contain the contribution from the 
quantum clusters of the model. Such exact results ob- 
tained in one-dimensional models containing both quan- 
tum and classical interactions can shed light on the prop- 
erties of the corresponding purely quantum ones and as a 
result into the more profound understanding of the ther- 
modynamic behavior of magnetic materials with a one- 
dimensional arrangement of exchange bonds. *^^^ For in- 
stance, certain qualitative properties of one-dimensional 
quantum spin systems can remain intact under the re- 
placement of some or even all quantum exchange inter- 
actions by Ising ones. In particular, if the magnetization 
curve has an intermediate plateau at some fixed value of 
the magnetization, the same feature may persist also in 
the modified system, although the region of parameters 
at which it appears can be different. ^'^^ In some cases, 
one even finds a rather good quantitative correspondence 
between the exact magnetization curves obtained for the 
modified system and that obtained by numerical methods 
for the underlying fully quantum system. ^^ 

Some real compounds may in fact be well approxi- 
mated by such mixed Heisenbcrg-Ising models. One ex- 
ample is given by a family of trimetallic coordination 
polymer compounds. ^^"^''^ Due to the presence of highly 
anisotropic Dy^'*" ions, these compounds yield a direct re- 
alization of low-dimensional spin systems with Ising and 
Heisenberg bonds. The single-chain magnet reported in 




FIG. 1. "Orthogonal-dimer" chain constructed from triangles 
of quantum s = 1/2 spins (large circles) interacting with each 
other via Heisenberg XX Z interactions (thick solid lines) and 
assembled to a chain by "classical" Ising spins (small circles). 
One coupling constant ( Ji) connects the quantum spins along 
a vertical line while a generally different coupling constant (J) 
connects the other two bonds in the quantum triangle. Fur- 
ther coupling constants are assigned to the Ising bonds con- 
necting left (Ki) and right (^^2) spins with the corresponding 
"classical" spins ai and CTi+i. 



Ref. 22 is an example of a Heisenbcrg-Ising chain of tri- 
angles with three-spin linear quantum clusters; another 
interesting molecular magnet (spin ring) based on the 
highly anisotropic properties of the Dy'^"'" ions has been 
reported in Ref. 23. Both systems have been analyzed 
by the transfer-matrix solution of Hciscnbcrg-Ising spin 
chains in Ref. 24. Another class of examples is given by 
chain magnets with two alternating magnetic ions, thus 
giving rise to alternating quantum and classical spins 
(see, e.g., Refs. 25-27). 

Here we continue to apply the exact transfer-matrix 
calculation for elucidating the thermodynamic and 
ground-state magnetic properties of one-dimensional spin 
systems with Ising and Heisenberg bonds"'^^"^^'^^ and per- 
form an analysis of the similarities and discrepancies be- 
tween magneto-thermal properties of the exactly solvable 
Heisenberg-Ising model and the underlying purely quan- 
tum Heisenberg system. We consider a one-dimensional 
s = 1/2 chain consisting of triangular clusters with XX Z 
Heisenberg interactions and single classical sites between 
them (see Fig. 1). Each spin at the single site is con- 
nected by bonds of Ising type with the two spins of its left 
triangle and with one spin of the right triangle. The cor- 
responding arrangement of sites into a chain is known as 
"orthogonal-dimer" chain. ^^"^^ However, since the sym- 
metry of the chain is already reduced by the nature of 
the interactions, here we also allow for more general in- 
teractions. 

The magnetization curve of the antiferromagnetic 
quantum orthogonal-dimer chain is famous for an infi- 
nite series of magnetization plateaux.'^" This infinite se- 
quence of plateaux is caused by the local conservation 
of the total spin on each vertical dimcr: for each fc, a 
ground state of the infinite system can be constructed by 
a periodic sequence of k vertical dimcrs in the triplet 
state s = 1 with two singlet dimers at the end. In 
the Heisenberg-Ising case considered here, only the main 
plateaux with magnetization equal to 0, 1/4, and 1/2 of 
the quantum orthogonal-dimer chain are observed in the 
zero-tcmpcrature magnetization curve. Depending on 
parameters of the models, we find various ground states 



with 4 and 8 sites per unit cell realized at these mag- 
netization values. Some of these ground states give rise 
to a doubling of the unit cell and a corresponding spa- 
tial modulation of the local magnetization. The appear- 
ance of such structures may be attributed to the left-right 
asymmetry for the triangular quantum clusters. Similar 
observations were previously reported in Refs. 16 and 17 
for the sawtooth chain with Ising and Heisenberg bonds. 
This paper is organized as follows. In section II we 
present the model and describe the exact computation of 
thermodynamic properties using a classical transfer ma- 
trix. In section HI wc then specialize to zero temperature 
and describe the ground states. Section IV contains a dis- 
cussion of the ground-state phase diagram in the presence 
of an external magnetic field. Section V presents some 
exact magnetization curves of the Heisenberg-Ising model 
for finite temperatures. Finally, in section VI, we present 
numerical results for the full quantum Heisenberg model 
and compare it to the Heisenberg-Ising system. In this 
comparison, we focus on two properties: the ground-state 
magnetization curve and the magnctocaloric effect. The 
paper ends with a conclusion. 



II. THE MODEL AND ITS EXACT SOLUTION 

We consider a spin chain consisting of s = 1/2 XX Z- 
Heisenberg triangular spin clusters, connected to each 
other via single spins situated between the triangles in 
such a way that each left-hand single spin is connected 
with the Ising bond CTrSf to the single spin of a triangu- 
lar cluster and each right-hand single spin is connected 
by the Ising bond to two spins in the vertical line of the 
triangular cluster (see Fig. 1). This chain is topologically 
equivalent to the famous orthogonal-dimer chain, ^^"'^^ al- 
beit with more general interactions. The Hamiltonian of 
the system is conveniently represented as a sum of the 
single-block Hamiltonians Hi 



H, 



^ = Er^'-^(^'+'^ 



i+l. 



(1) 



where 






V 
i,2 



Sis) 



+ Jl (Ai(S'j 2<5'i^3 + S^2Si,3) + 'S'i,2'S'i,3J 

-H2{Sli + SI2 + Sf 3) . 



(2) 



Here L is the number of blocks (cells) in the chain (the 
number of spins TV is 4L), Si^a stands for the spin opera- 
tors of the i-th triangle in the chain, a ~ 1,2,3 numbers 
the three sites of the triangle, and af is the z-component 
of the single spin situated at the left-hand side of the 
i-th triangle. Here, we allow the exchange interaction 
constants between the spins on the vertical line Ji, Ai 



and between the single spin of a triangle and spins on 
the vertical line J, A to be different as well as the inter- 
action constants for left and right single spins, A'l and 
K2, respectively. Magnetic fields acting on 5-spins and 
CT-spins are also considered to be different. 

We are going to derive exact expressions for all thermo- 
dynamic functions of the chain. To this end one should 
calculate the partition function 






(3) 



Spin variables of single spins can be regarded as classi- 
cal Ising ones, as they contribute only their diagonal z- 
component and hence can be replaced with their eigenval- 
ues ±1/2. Thus, in the partition function one should sum 
over all values of the ai and take a trace over the states 
of all spin operators S. Due to the special structure of 
interactions, the partition function can be evaluated ex- 
actly. To proceed, we first note that the Hamiltonians of 
different blocks commute 



[H„H,]=0. 



(4) 



Thus, one can expand the exponential in Eq. (3) and get 
the following expression 



(<t) *=1 



(5) 



where Tr.^ implies the trace over all states of the i-th 
block of the chain. The trace can be easily calculated 
once the block Hamiltonian Hi is diagonalized, which 
can be achieved analytically in our case. One finds the 
following eigenvalues 

Al,2 (^n CT^+l) == T (-^1 + 2J) 

±- {K,a, + 2K2(T,+i) T -H2, 
--.h{l + 2A^)±\Kia,T\H2, 



^^3,4 (ciiCTj+l) 
^5,6 



(o-j, (Ti+l) — - (JlAi - J - Q^a,,±a^+i] 

±-K2<Jt+i T ■^H2, 
±-K2a^+i T ■^H2, 



(6) 



where the following abbreviation is adopted: 

Qa,a' = v/8J2A2 + ( Ji(l - Ai) - J + 2{Kia + K2a')Y. 

(7) 
The corresponding eigenstates carry a definite value of 
S'^Q^ — SI + S2 + S^: the eigenvalues Ai^2 correspond to 
states with S^^^. ~ ±3/2 and the remaining eigenvalues 
to S^^^. = ±1/2. The eigenvectors read: 



Ai 
A2 

A3 

A4 

A5,7 



III 3 
' 2 

\V1 n 



A, 



6,8 



nt), 

^ (im>-int», 



x/2 



(iut>-un»: 



1 



(im> + int) + c±i;n», 



-2 



(u;t> + un>+5±im» 



(8) 

In each state | • • • ) the first, second, and third arrow de- 
notes an z-cigcnbasis of S'i.i, Si^2, and Si. 3, respectively. 
Furthermore, we have introduced the following abbrevi- 
ations 



c± 



c± 



Jl(l- 


-Ai)- 


-J+2{K2(TB^^ 


~ KiaL)±Q~ar,,an 


Jl(l- 


-Ai)- 


JA 
- J - 2{K2aB. - 


- KiaL)±Qar^-ari 



JA 



Here o'jk^l) stands for the value of the right (left) u-spin 
for the given triangle. 

After calculating the trace over the quantum degrees 
of freedom one arrives at the following expression for the 
partition function 



(<t) i=l 



O-JjCTi+i 



= TrT^ 



(9) 



where T is a 2 x 2 transfer matrix and Z^.^^._^-^ is the 
"partial" partition function of the i-th block: 

8 



-^"4^ cosh f ^3g2-A>.-2AV.^i 



Tr,e-'3«'=^e-'5^"("-'^'+^) 



Substituting A„ {ai,ai+i) from Eq. (6) one gets 

H2 - Kia, 



(10) 



Jl(l + 2Ai) 



cosh 



(11) 



+e' 



JlAi-J 



^f 2 cosh (3- 



'P 



^2 + ^2^i+l 



cosh I /3 



iai.-Oi+i 



According to Eq. (9), the partial partition functions for one block multiplied by the factor e'^^^'^'+'^'+i) yield the 



entries of the transfer matrix for a certain block: 



Z_ 



(12) 



Thus, calculation of the partition function for the system 
I 



under consideration boils down to evaluating the trace of 
the L-th power of a 2 by 2 transfer matrix in full analogy 
with the ordinary Ising chain (see for example Refs. 1 
and 32). The free energy per unit cell in the thermody- 
namic limit is obtained from the maximal eigenvalue of 
the transfer matrix (12): 



/ = In 




Z_ 



{{^ 



AZi 



Z_ 



(13) 



Now it is straightforward to obtain not only all thermody- 
namic properties, but also expressions for the sublattice 
and total magnetizations: 



tions read: 



«^-^Y.-^'s[T.'. 



-pn 



{<y) 



9HiJ B.H; 



Ms = 



3LZ 



2 / df 



(a) \i=l a=l / 



3 \dHy 



2//3.ffi 



M 



(cr) \i=l \ a=l J J 



\m.^\ms. 



(14) 



III. GROUND-STATE PROPERTIES 

From the eigenvectors for one block of the system one 
can construct all possible ground states of the whole 
chain. The ground states of the system under considera- 
tion can be classified by the values of magnetization, M, 
which we normalize to saturation values ±1. In partic- 
ular, one finds antiferromagnetic M ~ 0, ferrimagnetic 
(A/ ~ 1/2 or 1/4), and ferromagnetic M = 1 ground 
states. 

There are five "antiferromagnetic" ground states. 
Among them, there is a special one which is the closest 
analog of the fully dimerized ground state of the purely 
quantum orthogonal-dimer chain. ^^^'^'^ However, in con- 
trast to the latter, here the corresponding ground state 
exhibits a two-fold degeneracy per block due to the Ising- 
type interactions on the horizontal dimers. So, we will 
refer to this ground state as the degenerate antiferro- 
magnetic state, which we will abbreviate by |AF1). In 
the state |AF1), the spins on vertical dimers form a per- 
fect spin singlet, while the spins connected by a hori- 
zontal Ising bond are aligned antiparallel, but can freely 
change their overall orientation without changing the en- 
ergy of the ground state. The corresponding wave func- 



iAFi)=ni^±i,-T) 



(15) 



Here, ^ stands for the "down" ("up") orientation for 
the (7 spin from the i-th block, which appears simulta- 
neously with the \vi ^)i {\v_i ^)i) configuration of the 
ATATZ-triangle of the same block (compare Eq. (8) for the 
explicit form of Iwj-i a)i)- It is worth mentioning, how- 
ever, that the disentangled nature of the antiferromag- 
netic state of two spins connected by the Ising bond limits 
the analogy between the AFl and the perfectly dimerized 
ground state of the purely quantum orthogonal-dimer 
chain. In particular, the two-fold degeneracy per block 
in our case has its origin in the above mentioned differ- 
ence. Schematic figures illustrating the spin arrangement 
for the various ground states as well as their energies per 
block'^'^ can be found in Table I. 

There is another M = ground state with spatial pe- 
riod equal to the period of the system (four spins in the 
unit cell). This ground state corresponds to S^^^ = 1/2 
states of the quantum triangles which are aligned an- 
tiparallel with the s = 1/2 a-spins between them: 



|AF2)=ni 



.4). 



(16) 



Three further antiferromagnetic ground states are 
based on the eigenvectors of the triangular clusters which 
are asymmetric with respect to left and right a-spins. 
These ground states break the translational symmetry 
spontaneously, leading to a doubling of the unit cell 



L/2 



|AF3) = ni^'i,,-t)2.-i|t':i,,,i)2«, 

i=l 
L/2 

?.= 1 
L/2 

|AF5) = ni"l't)2.^iK'-|,;)2,. 



(17) 



TABLE I. Ground states of the system and their energies per block. '*'^ The function Q is defined in Eq. (7). Vertical ellipses 
denote the quantum singlet state of the corresponding spins; horizontal rectangles denote a disordered pair of two Ising spins 
with Stot = 0; gray triangles stand for the |i;']i ) ground states and a frame around a triangle indicates the degenerate 

superposition of \v^ ) and \v^ ). For more details compare the text. 
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Energy per block 



Spatial period 



Figure 
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undefined 
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The |AF3) ground state is characterized by an antiferro- 
magnetically ordered sublattice of triangles with S^^^ = 
±1/2 as well as by the same ordering of cr-spins. These 
two sublattices form an "up-up-down-down" overall ar- 
rangement of the local s = 1/2 magnetic moments. One 
can also notice that here we have an antiferromagnetic 
ordering of blocks with overall spin equal to ±1. In the 
|AF4) ground state, as in the previous case, one can see 
separate antiferromagnetic ordering of both cr-spins and 
triangles which are in the s = 1/2 state. However, inside 
each block they compensate each other (see Table I). In 
the |AF5) ground state one can see antiferromagnetic or- 
dering of fully polarized blocks. 

Now we turn to M ~ 1/4. Since this yields a frac- 
tional spin polarization per unit cell, Af = 1/4 ground 
states are accompanied by spontaneous breaking of trans- 
lational symmetry, ■^^'■^'* more precisely, doubling of the 
unit cell. In view of the spatial modulation of the lo- 
cal spin polarization, we will refer to the corresponding 
ground states as spin modulated ones: 



L/2 



\SMl)=Y[\vl^A)2^-l\vl^,i)2^, 

L/2 
|SM2)=ni"l,s't>2.-l|«i,a,;>2., 



(18) 



L/2 

\SM3)=Y[\vs,t)2^-l\vZ.,^,i)2^. 

These ground states have more complicated local spin 
arrangements than the antiferromagnetic ones. In con- 
trast to the latter, these "spin-modulated" ground states 
posses a doubled unit cell in which two blocks are not 
connected to each other by spin inversion (see Table I). 
Our model exhibits also three translationally invariant 
fcrrimagnetic ground states with M ~ 1/2: 



IFI) 



IF2) 



|F3) 



ni^i'^)* 



(19) 



ni 

L 

ni 

4 = 1 



ha>t)* 



A)^ 



Furthermore, there is also one family of ground states 
with a two-fold degeneracy per block with M = 1/2 in 
which all a spins point up, while the quantum triangles 
can freely and independently oscillate between \vi ^) and 

|w7 )■ However, these ground states appear only for the 

2 '^ 

highly symmetric choice of parameters A = Ai, J = Ji, 
Ki = K2- The corresponding eigenvectors are 



|F4) = n|yi,t). 



(20) 



TABLE II. Equations of the phase boundaries for J — Ji > 
0, Ki = 1x2 = K, and A = Ai (corresponding to the top 
panel of Fig. 2). The abbreviations q and q are introduced in 
Eq. (22). 



Phase 1 



Phase 2 



AF4 


SP 


h = 


AF4 


F4 


h = 


F4 


SP 


h = 


SM2 


F4 


h = 


AF3 


SM2 


h = 


SM2 


Fl 


h = 


AF3 


Fl 


h = 


Fl 


SP 


h = 


F4 


Fl 


K — 



Boundary 

T 



(4- A + 2K + q(-K)) 
i-3A + q{-K)) 
(2 -H A + k) 
(-3A-f 4K + g(K)) 
{-3A-2 + q{-n)) 
(8 4- Ai -4k + (j(-k)) 
(4- A-2Kg(-K)) 



1+|A 



where \Vi) can be either one of the states \vi g) or \vi 

2 2 ' 2 ' 

Finally, there is the fully polarized ground state 



isp)=ni^f't)4 



(21) 



IV. GROUND-STATE PHASE DIAGRAMS 

Ground-state phase diagrams can be derived from the 
energies of the ground states given in Table I. Since 
there are many free parameters, one needs to make some 
choices and it is natural to set certain exchange constants 
equal. We are going to focus on two cases. The first one 
corresponds to a situation with one exchange constant for 
all quantum bonds and another one for all Ising bonds 
(J = Ji = J, A'l = K2 = K) while the anisotropy con- 
stants for the two types of quantum bonds can still be 
different. The other case is well investigated in the con- 
text of the quantum orthogonal-dimer chain, ^^^-^^ where 
there are two interaction parameters, one for vertical and 
horizontal dimers, and another one for diagonal bonds 
{K\ = Ji = J,K2 = J = K). 



A. J ^Ji, Ki=K2=K 

Analyzing T = properties of the system under con- 
sideration one can find a rich variety of possible ground 
states. First of all, one can distinguish cases of ferro- 
magnetic and antiferromagnetic J. Hereafter, we will 
measure all coupling constants in units of \J\. 

In the antiferromagnetic region J > we observe a 
strong dependence on the ratio of A/Ai. For equal 
anisotropy in the triangles (A = Ai), the system exhibits 
an exceptional macroscopically degenerate ferrimagnetic 
ground state with M = 1/2, |F4). The top panel of Fig. 2 
presents the phase diagram for J > and A = Ai. This 
phase diagram has a point (K = 0, H = 0) where four 



-^ 



-i 







FIG. 2. Ground-state phase diagrams for antiferromagnetic 
Ji = J > and K\ = K2 = K. The top panel shows the 
case of isotropic interaction (A = Ai = 1). A characteristic 
feature is the appearance of the macroscopically degenerate 
ferrimagnetic phase F4 with Af = 1/2. The lower two panels 
show the case of Ising-like anisotropics, namely A = 0.8 > 
Ai = 0.5 (middle panel) and A = 0.5 < Ai = 0.8 (lowest 
panel). Note that A and Ai are the prefactors of the XY- 
term in the spin Hamiltonian. 



ground states merge, with a very high (eight-fold) degen- 
eracy per block. However, the origin of this degeneracy 
is trivial, as the system at such values of parameters is 



TABLE III. Equations of the phase boundaries for J = Ji > 
0, K\ = K2 = K, and A > Ai (corresponding to the middle 
panel of Fig. 2). 



Phase 1 



Phase 2 



AF4 

AF4 

F3 

AF3 

AF3 

AF3 

SM2 

SM2 

Fl 

F3 



SP 
F3 
SP 
F3 
Fl 
SM2 
F3 
Fl 
SP 
Fl 



Boundary 

h = 

h = 

h 

h 

h 

h 

h 



(4- Ai + 2K + g(-K)) 
(g(-Ac) - g(0)) 
(4 - Ai + 2ac + g(0)) 
(2k - q{0) + q{K)) 
{4-Ai-2k + q{K)) 
(-3Ai+g(K)) 
(3Ai+4K-2g(0) + g(K)) 
(8-h Ai -4K + g(K)) 



(Ai-g(O)) 



TABLE IV. Additional equations of the phase boundaries for 
J = Ji > 0, Ki — K2 = K, and Ai > A (corresponding 
to the lower panel of Fig. 2). The other equations of phase 
boundaries are the same as in Table III. 



Phase 1 



Phase 2 



AF4 

F2 

F2 

AFl 

SM2 

Fl 



F2 
SP 
AFl 
SM2 
F2 
F2 



Boundary 

ft = i (-3Ai + g(-K)) 
ft = I (2 -f Ai -f k) 

h=lK 

ft 
ft 



i(3Ai-g(K)) 

i (-3Ai + 4k -f g(K)) 



^l + |Ai 



just a set of triangles and single spins decoupled from 
each other. One also observes antiferromagnetic [M = 0) 
and ferrimagnetic spatially modulated [M = 1/2) ground 
states. The equations of the corresponding phase bound- 
aries are listed in Table IL Here and below we use the 
following notations: 



q{x) 

q{x) 
<1k{x) 
q;:{x) = v/8A2k2 + (A + 2xY 



V8A2 


+ (Ai + 2:e)2, 


V8A2 


+ (A + 2x)2 , 


v/8A2 


k2 + (Ai -1- 2a;)2 



(22) 



Throughout this section we use the abbreviations k ~ 
K/\J\ and/i = 7?/|J|. 

The middle panel of Fig. 2 presents the phase dia- 
gram for A/Ai > 1. Depending on the value of the 
ratio K = K/ J and the magnetic field measured in units 
of J [h = H/J), the system exhibits two antiferromag- 
netic, two ferrimagnetic, and one spin-modulated ground 
state, apart from the spin-polarized one. The present 
ground-state phase diagram differs from the previous one 
(see Fig. 2, top panel). Here, the degenerate ferrimag- 
netic ground state |F4) is replaced by |F3) and the region 
corresponding to |SM2) is smaller than in the previous 
case. An interesting feature of the symmetry breaking 
imposed by the distinct values of the anisotropy is the 
disappearance of the points where four phases merge. 
The equations for the corresponding phase boundaries 
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FIG. 3. Ground-state phase diagram for the case of ferro- 
magnetic Ji = J < 0, K\ = A'2 = /\, and A = 2, Ai = 3 
(XY-like anisotropy). 



TABLE V. Equations of phase boundaries for ferromagnetic 
J = Ji < and A-i = i<:2 = A" (Fig. 3). 

Phase 1 Phase 2 Boundary 

/i = i(Ai -4 4-2^ + 5(0)) 

h=\{2K + 2q{Q)-q{-K)-q{K)) 

h = \{2K~ 2q(0) + q{-K) + q(k)) 

/i = i (-6k + 2Ai - 8 + q(-K) + q{K)) 

h= i(Ai-4-2K-Fg(0)) 

fe = |k 



can be found in Table 111. 

The lower panel of Fig. 2 presents the phase diagram 
for the opposite relation between anisotropy constants 
A/Ai < 1 and antiferromagnetic J > 0. Here one new 
feature arises. Namely, one observes the macroscopi- 
cally degenerate antiferromagnetic ground state |AF1) 
between the |F2) and |SM2) states. Table IV contains 
the equations for all phase boundaries. 

Finally, for ferromagnetic J < one finds the phase 
diagram shown in Fig. 3. The corresponding equations 
of phase boundaries are listed in Table V. 

All phase transitions in Figs. 2 and 3 correspond to 
level crossings and thus should be considered as first- 
order transitions. 
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FIG. 4. Ground-state phase diagrams for the case J — K2 = 
K, J\ — A'l = J, which corresponds most closely to the fully 
quantum orthogonal-dimer chain. The top panel corresponds 
to isotropic interactions inside the triangles, A = Ai = 1, and 
antiferromagnetic J > 0. The middle (lower) panel shows the 
case of antiferromagnetic (ferromagnetic) J for the anisotropic 
case A = 2, Ai = 3 {XY-\ike anisotropy). 



The case where dimers have one coupling constant 
Ji = Ki = J and diagonal bonds another one J = K2 = 
K is the closest analog of the conventional purely quan- 
tum orthogonal-dimer chain. The ground state phase 
diagrams for this case are presented in Fig. 4. 

As before, one should distinguish ferromagnetic and 



antiferromagnetic regions of J. In both cases, the ap- 
pearance of different ground states is affected by the ra- 
tio A/Ai. The richest phase diagram exhibiting eight 
ground states corresponds to the isotropic case A = Ai 
and is shown in the top panel of Fig. 4. One observes 



TABLE VI. Equations of the phase boundaries for the case 
Ji = K\ = J > 0, J — A'2 = K and isotropic interaction 
inside triangles A — Ai = 1 (Fig. 4, top panel). 



Phase 1 Phase 2 



AF5 
AF5 
SM3 
AFl 
AFl 

F2 
AFl 
SM2 
AF3 
SM2 
SM2 

F3 
AF3 

F3 



SP 

SM3 
SP 

SM3 
F2 
SP 

SM2 
F2 

SM2 
F3 
F2 
F2 
F3 
SP 



Boundary 

u _ i 



i (1 + A + 2k) 
A + 2k) 
A + 2k) 



ft=i(3 

ft = I (2 

/i= |(-2-3A + 2AC + gr(«:)) 

ft = I (2 + 3A + 2k - 2Ag;7(0) + ^^(k)) 

h= |(2-3A + 2K + gr(fi:)) 

K.= 1 

/i=i(2K-Agr(0)+gr(«:)) 

h= 1(2- A + 4Ac + Agr(0)) 



A + 2k) 

3A-2K-gr(«:)) 
3A + 2K + gr(K)) 



TABLE VII. Equations of the phase boundaries for the case 
Ji = i^i = J > 0, J = 7^2 = if, and A = 2, Ai = 3 (Fig. 4 
middle panel). 

Phase 1 Phase 2 Boundary 



F3 

F3 

AFl 

F2 

AFl 

SM2 

SM2 

SM2 

AF3 

F3 

F2 



SP 

AFl 

F2 

SP 

SM2 

F2 

AF3 

F3 

F3 

SP 

F3 



h=l 

K = A 



(4k- Ai + 2 + g«(0)) 
(2 + 3Ai-q«(0)) 

(1 + Ai+2k) 

(2 + 3Ai -2K-gK(K)) 

(2-3Ai+2K + g«(K)) 

(-2-3A1 +2K + q«(K)) 

(2 + 3Ai + 2k - 2g«(0) + q^n)) 

(2k - g40) + q^K)) 

(l-Ai+5K + q«(0)) 

i/A 



TABLE VIII. Equations of phase boundaries for ferromag- 
netic Ji — Ki = J < and J = K2 = K (lower panel of Fig. 

4). 



Phase 1 



Phase 2 



F3 
AF3 
AF2 

AF2 



SP 

F3 

F3 

AF3 



Boundary 

|(4K + Ai+g40)) 



h 

/i = i (2k 

h 



i 



q.{0) + q4~K)) 
(2K-g«(0) + g„(-K-l)) 
(Ai - 1) 



and the antiferromagnetic |AF5) ground state is replaced 
by the ferrimagnetic state |F3). Still, there is a point 
of large entropy accumulation where four ground states 
|F3), |AF1), |F2), and |SP) become degenerate. In the 
symmetric but anisotropic case A = Ai 7^ 1, the value of 
the entropy at the corresponding point is S/L = 1.60944, 
corresponding to a five-fold degeneracy per block. Note 
also the line of degeneracy between the |F2) and |F3) 
states which one can see both in the top and middle panel 
of Fig. 4. This line is at K/J = Ai/A and corresponds 
to the degenerate |F4) ground state which in some sense 
is a superposition of |F2) and |F3). 

Lastly, the region of ferromagnetic J < is quite differ- 
ent. The lower panel of Fig. 4 demonstrates the ground- 
state phase diagram for ferromagnetic J and asymmetric 
anisotropic couplings A = 2, Ai = 3. The equations of 
the phase boundaries corresponding to the middle and 
lower panels of Fig. 4 are listed in Tables VII and VIII, 
respectively. 

Again, all phase transitions in Fig. 4 are first-order 
transitions in the sense that they correspond to level 
crossings. 



MAGNETIZATION PROCESS 



a broad region corresponding to the macroscopically de- 
generate antiferromagnetic ground state |AF1) which is 
the closest analog of the exact dimerized ground state 
of the purely quantum orthogonal-dimer chain. ^^^^° The 
system also exhibits four ground states with period dou- 
bling, |AF3), |AF5), |SM2), and |SM3). It is worth men- 
tioning that there is a point where the four ground states 
|SM3), |AF1), |F2), and |SP) become degenerate, which 
implies a large entropy accumulation. For the isotropic 
case (A = Ai = 1) presented in the top panel of Fig. 4, 
the position of this point is K/J = —1/2, H/J = 1/2 and 
the corresponding entropy per unit cell is S/L — 1.44364. 
The equations of the phase boundaries can be found in 
Table VI. Throughout this section we use the abbrevia- 
tions K = K/\J\ and h = H/\J\. 

The structure of the phase diagram for anisotropic in- 
teractions and antiferromagnetic J, K > is qualita- 
tively similar. The middle panel of Fig. 4 shows the 
ground state phase diagram for A = 2, Ai = 3, and J > 
0. Here one can see that, in comparison with the isotropic 
case, there is no region corresponding to the |SM3) state 



The phase diagrams show that the model (1), (2) ex- 
hibits a large diversity of magnetic behaviors. Namely, 
depending on the values of parameters, the following 
step-like transitions between plateaux in the magnetiza- 
tion curves can be observed when the magnetic field is 
varied from zero to the saturation value: M = — > 1, 
A/ = ^ 1/2 -^ I, M = ^ 1/4 -^ 1/2 ^ 1, A/ = 
0^1/4^ 1, AI = 1/2 -^ 1, and M = 1/4 -^ 1/2 -^ 1. 
As there are no excitation bands in the model, all transi- 
tions between different ground states are strictly step-like 
at T = 0. 

Here we illustrate some variants of the magnetization 
process and demonstrate the effect of a finite temper- 
ature using the exact solution of the model. For suffi- 
ciently low temperatures, the form of the magnetization 
curve almost coincides with the zero-temperature limit. 
Further details for T ~ will be given in the comparison 
with the full quantum mechanical Heisenberg model in 
Section VI B below. 

Fig. 5 shows the magnetization curve for the most sym- 
metric case J = Ji = Ki = K2 = 1 and A = Ai = 1, 
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FIG. 5. Magnetization curve for Ji = Ki — 1, J — K2 = 1, 
and A = Ai = 1. The corresponding ground states can be 
seen both in the top panel of Fig. 2 and the top panel of Fig. 
4. For very low temperature (solid line) the sequence of step- 
like transitions between the corresponding ground states is 
recovered. Dashed and dotted lines correspond to T/J — 0.05 
and T/J = 0.1, respectively. 



FIG. 7. Magnetization curve for Ji — K\ = 1, J = K2 ~ 
1.19034, and A = 2, Ai = 3 {XY-like anisotropy). The solid 
line corresponds to T/J = 0.001 and shows sharp transitions 
between three plateaux: M = 1/4, 1/2, and 1. Dashed and 
dotted lines correspond to T/J = 0.05 and T/J = 0.1, re- 
spectively. 
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only two intermediate plateaux at M = and 1/2, re- 
spectively. For the values of parameters corresponding 
to the figure (J = Ji = 1, Ki = K2 = —1 and A = 0.5, 
Ai = 0.8), these plateau correspond to the |AF4) and 
|F2) ground states (see lower panel of Fig. 2). 

Fig. 7 shows the magnetization curves for Ji = Ki = 1, 
J = K2 = 1.19034 and A = 2, Ai = 3 which corresponds 
to the middle panel of Fig. 4. This exceptional value of 
K/ J corresponds to the point of the phase diagram where 
the three ground states |AF1), |SM2), and |F3) become 
degenerate, {K/J = t^( 



= 0) and the 
first plateau in the magnetization curve corresponds to 

M = 1/4. 



VI. HEISENBERG MODEL 



FIG. 6. Magnetization curve for J = Ji — 1, K\ = K2 ~ —1, 
and A = 0.5, Ai — 0.8 (Ising-like anisotropy). The solid line 
corresponding to T/J — 0.001 shows sharp transitions be- 
tween three plateaux: M = 0, 1/2, and 1. Dashed and dotted 
lines correspond to T/J = 0.05 and T/J = 0.1, respectively. 



for T/J = 0.001, 0.05, and 0.1. One can see here a 
sequence of magnetization plateaux corresponding to the 
following transitions between the ground states at T = 0: 
|AF3) -^ |SM2) -^ |F4) ^ |SP) which corresponds to 
both phase diagrams from Fig. 2 and Fig. 4. The line 
K/J = 1 from the upper panel of Fig. 4 where |F2) and 
|F3) become degenerate corresponds to the degenerate 
|F4) ground state. 

Fig. 6 exhibits another form of magnetic behavior with 



In the following we would like to test to which extent 
the reduction to Ising interactions for certain bonds af- 
fects the properties of the model. For this purpose we 
will use the Heisenberg-Ising variant of the model (1), 
(2), i.e., we will put A = Ai = 1. 

The comparison will be carried out with the Heisenberg 
analog of the Hamiltonian 



N/i 



i.2 -t- Ji.i) 



Jl Si,2 ■ S 



i,3 



+Ki Si,i ■ Si^4 + K2 {Si,2 + S'l.a) • S'i+1^4 
-hJ^SiX (23) 

a=l ' 

Note that we have replaced the Ising spin oi in Eqs. (1), 
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(2) with a full s = 1/2 Heisenberg spin Si,4- 



A. Method 

Firstly, one can use conservation of the total 5*^ = 
Tliia^ia to relate the energy E{S^,H) in a magnetic 
field to that at zero field 



E{S\H)^E{S'',0)-HS'' 



(24) 



It is therefore sufScient to diagonalize (23) in a given 
sector of 5^ for iJ = 0. 

Secondly, we can introduce composite operators 



Ti — Si^2 + Si^z 
and rewrite the H ~ case of (23) as follows: 



(25) 



N/4 



1 • 



Ji Tf-- 



+Ki S.,1 • S^A + A'2 T, ■ S,+i,4 1 . (26) 

This puts our model into the class of one-dimensional 
models with local conservation laws.-^^'^^'^®'^®"'^^''^^"^^ 
This property can be used to simplify the diagonaliza- 
tion: Since the total spin Tj on each vertical dimer is 
separately conserved, diagonalization of (23) is reduced 
to diagonalization of (26) for all possible combinations 
T; = 0, 1 (with Tf = T,{Ti + l)). In fact, a singlet 
Ti ~ splits the chain into "fragments": all states for 
a given singlet pattern are product states of the corre- 
sponding eigenstates of the open fragments separated by 
these singlets. Hence, it is sufficient to diagonalize all 
open fragments with up to N/A— 1 blocks in which Ti = 1 
and one periodic system with N/A with all T^ = 1 in or- 
der to obtain the spectrum for a periodic chain with N 
sites. ^^'^^ For K2 = J, the use of reflection symmetry al- 
lows us to fully diagonalize a system with N = 24 spins 
1/2, a task which in view of the large unit cell of the 
model would be very demanding if not impossible for the 
original Hamiltonian (23). 

Ground states consist of a periodic repetition of one 
fragment of a certain length or a periodic system of a 
given length with all Ti = 1.'^° In the former case, we can 
construct the exact thermodynamic limit of the magneti- 
zation curves, provided that we have access to sufficiently 
long fragments. For the periodic system, we use (i) val- 
ues for plateau boundaries obtained at sufficiently large 
N and (ii) the standard midpoint method^®'^^ as applied 
to the largest available system size to approximate the 
thermodynamic limit in smooth regions of the magneti- 
zation curve. 

The main limitation for the computation of ground- 
state properties is the diagonalization of a periodic sys- 
tem with all Ti = 1. To push this a bit further, we 
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FIG. 8. (Color online) Zero-temperature magnetization 
curves for two cases where all exchange constants are an- 
tiferromagnetic: Ji — K\ = 1, J = K2 ~ 0.7 (a) and 
Ji = J — K\ = K2 = 1 (b). Dashed lines show numeri- 
cal results for a Heisenberg chain with A'^ = 40 spins 1/2 and 
thick full lines an extrapolation to the thermodynamic limit. 
Thin full lines show the exact solution of the Heisenberg-lsing 
model in the thermodynamic limit A*' = 00. The inset of panel 
(b) shows the finite-size behavior of the edge of the M = 
plateau (spin gap) and the two edges of the M = 1/2 plateau 



have used the ALPS 1.3 implementation^'' of the density- 
matrix renormalization group (DMRG) method. ^^^^^ Be- 
low we show results based on DMRG for A^ = 40 in the 
regime M < 1/2 and iV = 64 with m = 400 kept states 
as well as A^ = 96 with m = 600 kept states. In the case 
Ji = J = 1, A'l = A'2 = — 1, we have exceptionally in- 
creased the number of kept states to m = 600 and 800 for 
A^ = 64 and 96, respectively. We have always also per- 
formed computations with lower m to ensure convergence 
of the data. 
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TABLE IX. Local spin polarizations on the M — 1/4 plateau 
for Ji = /<i = 1, J = A'a = 0.7 (see Fig. 8(a)). 





Heisenberg 


Heisenberg-Ising 


Sli 


0.3312 


0.5 


s!,2 + si3i=Tn 








Si+IA 1 '^'+1 


0.3312 


-0.5 


•Ji + l,! 


-0.2237 


-0.1022 


'-'i'+l,2 + '-'i+1,3 (= 2^ + l) 


0.7849 


0.6022 


'S'i'+2,4 / ^^^+2 


-0.2237 


0.5 



B. Magnetization curves 

Now we present examples of numerical results for mag- 
netization curves of the Heisenberg model (23) and com- 
pare them with the Heisenberg-Ising model (1), (2) with 
A = Ai = 1. 



1. ,h ^Ki = l, J^K2= 0.7 

Fig. 8(a) shows the case Ji = A'l = 1, J = A'2 = 
0.7 which has been discussed prcviously.^^^'^^ Our new 
exact diagonalization results agree with the previous 
ones. Note that the spin-1/2 Heisenberg model most 
likely gives rise to an infinite sequence of magnetization 
plateaux which accumulate just below M = 1/2.'^" This 
infinite sequence of plateaux is not reproduced by the 
Heisenberg-Ising model (red thin full line in Fig. 8(a)) 
which does, however, reproduce some of the most promi- 
nent features: (i) a spin gap at M = 0, (ii) two broad 
plateaux at M = 1/4 and 1/2, and (iii) a jump from 
M = to 1/4. This choice of parameters corresponds to 
a cut through the top panel of Fig. 4 at K/J = 0.7 such 
that one can read off the corresponding ground states of 
the Heisenberg-Ising model. 

For M = 0, the ground state of the Heisenberg model 
is characterized by all Ti = 0. This leaves isolated hori- 
zontal dimers which also form singlets. The correspond- 
ing ground state of the Heisenberg-Ising model is the 
|AF1) state. This is in fact the closest analog of the 
fully dimerized ground state: the vertical dimers are also 
in the singlet state for the |AF1) state. However, due to 
the missing quantum fluctuations the horizontal dimer 
bond is unable to form a singlet. The two local classical 
antiparallel spin configurations appear both with equal 
probability in the |AF1) state instead of forming a sin- 
glet superposition. This also gives rise to a macroscopic 
ground-state degeneracy of the Heisenberg-Ising model 
with entropy S/L = In 2 per unit cell. 

For M = 1/4, Ti = and 1 alternate in the ground 
state of the Heisenberg model, i.e., the ground state 
breaks translational symmetry spontaneously. The same 
translational symmetry breaking is observed in the cor- 
responding |SM2) state of the Heisenberg-Ising model. A 
more detailed comparison of the two states is given by the 
local spin polarizations along the field direction shown in 



Table IX. For the Heisenberg model, this follows from the 
local expectation values in the S^ = 1 sector of a frag- 
ment consisting of a single block. For the Heisenberg- 
Ising model one needs instead the expectation values in 
the |w7 ) state of a single quantum triangle. Using the 

explicit form of the eigenvector given in Eq. (8) one ar- 
rives at 



KJ(^2+^I)K, 



{vljstlvj 



(27) 






Combination with the |SM2) wave function given in 
Eq. (18) yields the corresponding entries in Table IX. 
The entries in Table IX exhibit a similar structure of the 
M = 1/4 state at Ji = Ki = 1, J = K2 = 0.7 in the 
Heisenberg and Heisenberg-Ising models. In particular, 
vertical dimers in the singlet state alternate with pre- 
dominantly polarized ones. The quantitative differences 
observed in Table IX can be attributed to the interac- 
tion of the Heisenberg-Ising model breaking the spatial 
reflection symmetry of the Heisenberg model. 

In the AI = 1/2 state of the Heisenberg model one 
flnds all vertical dimers in the state Ti = 1. This state is 
non-degenerate and thus translationally invariant. The 
spin polarizations in this state have to be determined nu- 
merically; using DMRG one finds Si-^ 



5^4 = 0.0640 



and SI2 + SI3 = TI ^ 0.8721.^^ The corresponding 
ground state of the Heisenberg-Ising model is the |F2) 
state. Accordingly, the spin polarizations are read off 
from Eq. (19) as Sl^ = Oi = 1/2 and 5f^2 + ^Iz = 0- So, 
in this case the difference between the Heisenberg model 
and the Heisenberg-Ising model is bigger although the 
AI = 1/2 ground states still have the same translational 
symmetry. 



2. .h = J = Ki=K2 = l 

The effect of increasing J = K2 from 0.7 to 1 is 
demonstrated in Fig. 8(b). The magnetization curve 
of the Heisenberg-Ising-model (which was in fact al- 
ready presented in Fig. 5) remains similar to the one 
in Fig. 8(a). However, some of the plateau states of 
the Heisenberg-Ising-model change; the sequence is now 
|AF3) -> |SM2) -> |F3) -> |SP) (see top panel of Fig. 4). 
The Heisenberg model still has a spin gap at M = 0, but 
its state is also changed: now it arises out of a periodic 
system with all T^ = 1, as already observed in Rcf. 44. 
The M = 1/2 state of the Heisenberg model also corre- 
sponds to all Ti = I and is non-degenerate. Since in this 
case finite-size effects are important at M = 1/2, we have 
performed a finite-size analysis which is shown in the in- 
set of Fig. 8(b). While the width of the M = 1/2 plateau 
initially shrinks with growing A^, A^ = 96 is a good ap- 
proximation to the thermodynamic limit 1/A^ = 0. In 
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TABLE X. Values of local spin polarizations for the ground states corresponding to the magnetization plateaux in Fig. 9. The 
|AF4) ground state breaks translational symmetry such that blocks with different signs alternate. Hence, the upper (lower) 
sign should be taken for an even (odd) cell of the |AF4) state. 



Case 


Heisenberg-Ising 






Heisenbcrg 






Ground-state label Sf i 


0"i 


Si, 2 + •S'j's 


Sli 


'S'i',4 


Tf 


= 5f,2+Sf,3 


Figs. 9(a), (b), M = 


AF2 -0.4621 


1/2 


0.9621 


— 


- 




— 


Fig. 9(c), M = 


AF4 =Fl/6 


±1/2 


±2/3 


- 


- 




- 


Fig. 9(a), M = 1/2 


F3 -1/6 


1/2 


2/3 


0.2640 


0.2640 




0.4720*^^ 


Fig. 9(b), Af = 1/2 


F2 1/2 


1/2 





1/2 


1/2 







Fig. 9(c), M = 1/2 


F4 1/6 


1/2 


1/3 


-0.0799 


0.4287 




0.6512'^^ 



particular, we conclude that the M — 1/2 plateau sur- 
vives for N = oo. All further plateaux in the magnetiza- 
tion curve of the Heisenberg model are probably gone. 



3. A'l < 

Fig. 9 shows three examples of magnetization curves 
with ferromagnetic Ki < and potentially a further fer- 
romagnetic exchange constant. In all three cases shown 
in Fig. 9, the Heisenberg-Ising model has plateaux at 
M = and 1/2, separated by jumps in the magnetiza- 
tion curve. These two plateaux are also present in the 
Heisenberg model. However, in the case Ji = J = 1 and 
Ki = K2 = — 1 a finite-size analysis is again necessary 
to show that the spin gap (corresponding to the AI = 
plateau) persists in the thermodynamic limit. According 
to the data shown in the inset of Fig. 9(c), we expect 
A^ = 96 to be a good approximation to the thermody- 
namic limit even if the spin gap decreases substantially 
with growing A^ for the smallest values of A^ considered. 

At a quantitative level, the spin gap of the Heisen- 
berg model compares favorably to that of the Heisenberg- 
Ising model in the case of Fig. 9(a) and the width of 
the M = 1/2 plateau compares favorably in all three 
cases presented in Fig. 9. The best agreement is ob- 
served for the case presented in Fig. 9(b) (Ji = J = 1, 
Ki = K2 = —1). In this case, the jump between 
M = 1/2 and 1 is reproduced exactly and both mod- 
els give rise to a jump at the lower edge of the M = 1/2 
plateau, albeit not exactly at the same magnetic field. 

For the three examples shown in Fig. 9, the spin gap 
of the Heisenberg model always arises out of a periodic 
pattern Ti = 1. This conclusion agrees with Ref. 44 in 
the case of Figs. 9(a) and (b). In fact, the behavior at 
low magnetic fields is identical in these two cases, includ- 
ing the feature at low magnetizations in Fig. 9(b). To the 
best of our knowledge, the case presented in Fig. 9(c) has 
not been discussed previously. The ground state of the 
Heisenberg-Ising system at M = is the |AF2) one in 
the case of Figs. 9(a) and (b) and the |AF4) one in the 
case of Fig. 9(c). The corresponding local spin expecta- 
tion values are given in Table X. In the Heisenberg model 
at zero magnetization, individual spin expectation values 
vanish for symmetry reasons. Hence, for the Heisenberg 



model at M = 0, we have to look at correlation func- 
tions in order to determine the structure of the ground 
state. Fig. 10 shows selected correlation functions com- 
puted by exact diagonalization for a periodic A^ = 32 
system as a function of distance j between cells. The 
correlations of the left panel reveal a structure which is 
similar to the |AF2) state, i.e., the spins S^^ and S^^ 
next to the composite spin T^ point in the opposite direc- 
tion of this composite spin while spins of the same type 
show a tendency of aligning parallel between different 
cells. By contrast, in the right panel of Fig. 10 one ob- 
serves a tendency of antiparallel alignment between cells 
at an odd distance while the internal structure of a cell 
is the same as before, as expected for the |AF4) state. 
In the case of the left panel of Fig. 10, correlations are 
short-ranged while for the right panel they decay more 
slowly in accordance with the size of the spin gap (given 
by the width of the M = plateau) in Figs. 9(a), (b) and 
Fig. 9(c), respectively. Overall, we conclude that for all 
three cases considered in Fig. 9, the M = ground state 
of the Heisenberg-Ising model is a good qualitative rep- 
resentation of the M = ground state of the Heisenberg 
model. 

Turning to finite magnetizations, we observe that in 
the case of Figs. 9(a) and (c), the M = 1/2 plateau state 
of the Heisenberg model arises out of a system with all 
Ti = I while in the case of Fig. 9(b) it is characterized 
by all Ti = 0. The corresponding local spin expectation 
values are given in Table X. The ground states of the 
Heisenberg-Ising model corresponding to the AI = 1/2 
plateau in the three cases presented in Figs. 9(a), (b), 
and (c) are |F3), |F2), and |F4), respectively. The spin 
expectation values for these states are included in Ta- 
ble X. Note that the |F4) state is macroscopically de- 
generate (see Eq. (20)) such that the average over this 
manifold yields effectively the average of the local expec- 
tation values in the |F2) and |F3) states. In the case of 
Fig. 9(a), the reflection-symmetry breaking interactions 
of the Heisenberg-Ising model lead to a breaking of reflec- 
tion symmetry in the M = 1/2 state which is not present 
in the Heisenberg model. Otherwise, the structure of the 
M ~ 1/2 states of the Heisenberg-Ising and Heisenberg 
model is at least similar in the case of Figs. 9(a) and (c). 
In the case of Fig. 9(b) we even find an identical M =1/2 
state as reflected by identical local spin expectation val- 
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FIG. 9. (Color online) Same as Fig. 8, but for cases which 
contain ferromagnetic exchange constants, namely Ji = Ki = 
-1, J = 7^2 = 1 (a), Ji = J = K2 = I, Ki = -I (b), and 
Ji = J = 1, Ki = K2 = —1 (c). The inset of panel (c) shows 
the finite-size behavior of the spin gap (edge of the M — 
plateau) at Ji = J = 1, A'l = K2 = — 1- 
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FIG. 10. (Color online) Correlation functions of the z- 
components for the Heisenberg model with A'' = 32 sites 
at M = 0. j measures the distance in units of cells; due 
to the periodic boundary conditions, the maximum accessi- 
ble distance is A''/8 = 4. The left panel is for Ji = ±1, 
J = K2 = 1, Ki = — 1 whereas the right panel shows results 
for Ji = J ^ 1, Ki ^ K2 ^ -1. 



C. Entropy and cooling rate 

We now turn to thermodynamic properties. Since one 
of the most important properties of the present type of 
models is the magnetocaloric effect, ^'^'^^"^^'^^'^^"^^ we fo- 
cus on the entropy S* as a function of magnetic field H 
and temperature T as well as the adiabatic cooling rate 
as a function of magnetic field. These quantities have 
been computed for the Heisenberg model by full diago- 
nalization. 

The orthogonal arrangement of dimers in our model is 
very similar to the two-dimensional Shastry- Sutherland 
model. ^^ The latter model is found to have a good re- 
alization in SrCu2 (603)2 with a ratio between dimer 
and plaquette couplings estimated to be close to 0.7 (see 
Rcf. 70 for a review). Hence, the case Ji = Ki = 1, 
J ~ K2 = 0.7 can be viewed as the one-dimensional coun- 
terpart of the two-dimensional Shastry-Sutherland model 
for SrCu2(B03)2- Indeed, in this parameter regime, both 
the one- and two-dimensional Heisenberg model have an 
exact dimer ground state (compare section VI B and 
Ref. 69, respectively) and it was noted early on that 
the magnetization curve of SrCu2 (603)2 increases very 
steeply above a magnetic field H w 20T, like a system 
of isolated dimers^^ or the one-dimensional model, see 
Fig. 8(a). 

Fig. 11 shows our results for the entropy per spin 
S/N a.t Ji = Ki = 1, J = K2 = 0.7 as a function 
of magnetic field H and temperature T. Here we fo- 
cus on magnetic fields corresponding to magnetizations 
M < 1/4. At a quantitative level there are certain dif- 
ferences between the Heisenberg-lsing model (top panel 
of Fig. 11) and the Heisenberg model (bottom panel). 
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FIG. 11. (Color online) Entropy per spin S/N as a function of 
H and T for Ji = A'l = 1, J = A'2 = 0.7. The top panel shows 
results for the Heisenberg-Ising model in the thermodynamic 
limit and the lower panel the Heisenberg model with N — 24 
spins. Lines show values of constant entropy, starting with 
S/N = 0.02 at the bottom and increasing in steps of 0.02. 



FIG. 12. (Color online) Entropy per spin S/N as a function of 
H and T for Ji — J — K2 ~ 1, Ki = — 1. The top panel shows 
results for the Heisenberg-Ising model in the thermodynamic 
limit and the lower panel the Heisenberg model with A'" = 24 
spins. Lines show values of constant entropy, starting with 
S/N = 0.05 at the bottom and increasing in steps of 0.05. 



Firstly, from the magnetization curve Fig. 8(a) we al- 
ready know that the transition field He between the 
M = and 1/4 states is shifted. Secondly, we also 
know that the M = state of the Heisenberg model is 
non-degenerate whereas the Heisenberg-Ising model has 
a macroscopic degeneracy giving rise to a residual en- 
tropy S/N — In 2/4 = 0.173..., and accordingly, the 
values of S differ in the low-field low-temperature regime 
of Fig. 11. Nevertheless, the Heisenberg-Ising model and 
the Heisenberg model share important qualitative fea- 
tures. In particular, in both cases another local exci- 
tation comes down as H — )■ He- Noting that in the 
Heisenberg model these local excitations are not allowed 
to occupy two consecutive blocks for M < 1/2, one can 
count the degeneracy at H = He using a 2 x 2 trans- 
fer matrixi'32.64 and finds S/N = ln((l + \/5)/2)/4 = 
0.120 ... in the limit T ^ 0. A similar 3x3 transfer- 
matrix procedure for the Heisenberg-Ising model yields 
S/N = ln(l + ^/2)/4 == 0.220... = In 2/4 + 0.047. . . 
for the zero-temperature limit at H = He- Recall also 
that the ground state at M = 1/4 is two- fold degen- 
erate. For the A^ = 24 system shown in the lower 
panel of Fig. 11, this gives rise to a finite-size value 
S/N = In 2/24 = 0.028 . . . ior H > He and in the limit 

The additional residual entropy at H — He gives rise to 



substantial cooling during adiabatic (de)magnetization, 
as also demonstrated by the constant entropy curves 
in Fig. 11. In fact, if one starts with a sufficiently 
low temperature at H = 0, one reaches T = as 
H approaches He, both in the Heisenberg model and 
the Heisenberg-Ising model. Similar behavior is sug- 
gested in the two-dimensional Shastry- Sutherland model 
by the flat triplet branch above the dimerized ground 
state. ^^ Indeed, there is indirect evidence for substantial 
cooling in SrCu2 (603)2 during pulse-field magnetization 
experiments.^'^ 

The other case of special interest is Ji = J = K2 = 1, 
Ki = —1. Not only does Fig. 9(b) demonstrate quanti- 
tative agreement between the Heisenberg model and the 
Heisenberg-Ising model at the saturation field, but this 
case can also be considered as a realization of the hard- 
monomer universality class of localized magnons.^^"'*^ 

Fig. 12 shows the entropy per site S/N for the 
Heisenberg-Ising model (top panel) and the Heisenberg 
model (lower panel) at Ji = J = K2 = 1, Ki = —1. 
The global behavior of both models is similar. The main 
qualitative difference can be observed at small magnetic 
fields when the zero-field gap is closing: the Heisenberg- 
Ising model has a macroscopic ground-state degeneracy 
S/N = In 2/4 at the position of the first step in the 
magnetization curve Fig. 9(b) whereas in the Heisenberg 
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FIG. 13. (Color online) Adiabatic cooling rate as a function 
of magnetic field H for three different temperatures at Ji = 
J = K2 = 1, K\ = —1. Results for the Heisenberg-lsing 
model and hard monomers are for the thermodynamic limit 
while those for the Heisenberg model are for N = 24: spins. 



model the degeneracy is lifted in correspondence with the 
smoother transition. Nevertheless, also the Heisenberg 
model exhibits an enhanced low-temperature entropy as 
the r = magnetization increases from M = to 1/2. 
This enhanced entropy is reflected by the minimum in 
the constant entropy curves in the lower panel of Fig. 12. 

Close to the saturation field H = 2 wc find even 
quantitative agreement between the Heisenberg-lsing and 
Heisenberg model. In particular, both models now give 
rise to a zero-temperature entropy S/N = In 2/4 at 
H ^ 2 such that one could cool to arbitrarily low tem- 
peratures by adiabatic magnetization or demagnetiza- 
tion when approaching the saturation field from below 
or above, respectively. Note that in this regime finite- 
size effects are very small, as is evident, e.g., from the 
effective hard-monomer description. ^^"^^ Hence, we can 
use our N = 2A exact diagonalization data for reference. 

For a more detailed comparison we use the adiabatic 
cooling rate which can be written as 



dr 

dH 



dS/dHJT 
dS/dT\H 



(28) 



This is nothing else but the slope of the constant entropy 
curves in Figs. 11 and 12. The adiabatic cooling rate 
has the advantage over the entropy that it is directly 
accessible in experiments. ''^"^^ 

For hard monomers, the right-hand side of Eq. (28) 
can be evaluated easily using the free energy per unit cell 
/ = -T In (1 + cxp{-{H - H^)/T)) and S = -d f/dT: 



dT 
dH 



monomer 



T 



H - Hr 



(29) 



Fig. 13 compares the cooling rate for the Heisenberg 
model, the Heisenberg-lsing model, and the effective hard 



monomer low-energy description for the same parameters 
as in Fig. 12 for three different temperatures. All three 
descriptions exhibit a strongly enhanced cooling rate for 
H Ri He = 2 which reflects the T = entropy at He = 2 
(see Fig. 12). Not surprisingly, the quantitative agree- 
ment between the different descriptions is best at the 
lowest temperature T = 0.01. The main difference be- 
tween the Heisenberg-lsing model and hard monomers is 
that the singularity in Eq. (29) at H = He and any T > 
is removed in the Heisenberg-lsing model by the presence 
of further degrees of freedom at higher energies. This im- 
proves the agreement with the full Heisenberg model at 
higher temperatures. 



VII. CONCLUSION 

We have presented a detailed analysis of the ground 
states and magnetic properties of a one-dimensional lat- 
tice spin model with a structure which is very similar 
to the famous orthogonal dimer chain, ^^"'^^ although the 
exchange constants arc more general. In the Heisenberg- 
lsing variant, the main ingredient is the block structure 
of the Hamiltonian. Each block consists of a triangle of 
s = 1/2 spins interacting with each other via an XXZ- 
interaction and of one single spin, which is connected to 
two "dimer" spins from the triangle by Ising bonds. This 
particular structure renders the system exactly solvable 
by the classical transfer-matrix method. 

First, we investigated the large variety of T = ground 
states in an external magnetic field. Some of the ground 
states break translational symmetry spontaneously, thus 
giving rise to doubling of the unit cell. Moreover, 
we found two macroscopically degenerate ground states. 
One of them, the zero-magnetization state |AF1), is the 
closest analog of the famous dimerized ground state of the 
underlying purely quantum orthogonal dimer chain. ^^^'^'^ 
The main difference is that the Heisenberg-lsing system 
cannot form a quantum dimer on the horizontal bond, 
but instead the two classical antiparallel spin configura- 
tions are equally likely to be found on each given horizon- 
tal bond, thus yielding a two-fold degeneracy per block. 
The other ground state with a two-fold degeneracy per 
block, |F4), carries magnetization M = 1/2 and has a dif- 
ferent origin. It appears only at very symmetric values of 
parameters and corresponds to a degeneracy between the 
symmetric and antisymmetric S^^^ = 1/2 ground states 
for the two coupled quantum spins on a vertical dimer. 

The second focus of our investigation was a compari- 
son of the Heisenberg-lsing model with numerical results 
for the full Heisenberg model. Certain features of the 
full Heisenberg model are not present in the Heisenberg- 
lsing variant. For example, for certain antifcrromagnetic 
values of the exchange constants, the Heisenberg model 
also exhibits ground states with a periodicity larger than 



two 



30,44 



and rclatedly an infinite sequence of plateaux 
in its magnetization curve. '^'^ Only the main plateaux at 
Af = 0, 1/2, and 1/4 are also present in the Heisenberg- 
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Ising variant. 

When considered as an approximation to the Heisen- 
bcrg model, the Heisenberg-Ising model performs in gen- 
eral better if some of the Ising exchanges are ferromag- 
netic. In the three examples with a ferromagnetic A'l < 
which we have discussed, the presence and nature of 
plateaux at M = and 1/2 is qualitatively reproduced 
to the extent possible. Between these plateaux, the 
Heisenberg-Ising model always yields jumps in the T = 
magnetization curve while the Heisenberg model may ex- 
hibit smooth transitions. Also the critical fields can differ 
noticeably. Remarkably, in the case Ji = J = K2 = 1, 
Ki = —1 we found that the Heisenberg-Ising model re- 
produces results for the Heisenberg model even at a quan- 
titative level close to the saturation field (see Figs. 9(b) 
and 13). This exact correspondence may be attributed 
to the product-state structure of the low-energy states of 
the Heisenberg model in this parameter regime, *^^^^ and 
we speculate that it is generic for a sufficiently strong 
coupling of the vertical dimer exchange Ji . 

We believe that further continuation of recent devel- 
opments in the field of magnctochcmistry^^"^'^ will make 
it possible to obtain magnetic materials which match 
the model considered in the present paper. At the mo- 



ment, the most direct application of our system probably 
is as a one-dimensional version of the two-dimensional 
Shastry-Sutherland model for SrCu2(B03)2.^° In par- 
ticular, we have observed an enhanced magnetocaloric 
effect upon closing the zero-field spin gap, both in 
the one-dimensional Heisenberg as well as its simpli- 
fied Heisenberg-Ising variant, in accordance with indi- 
rect evidence for cooling by adiabatic magnetization of 
SrCu2(B03)2.'' 
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